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ABSTRACT 

We investigate the heat source of the neutral gas comprising damped Lya systems. Unlike the Lya 
forest, where the extragalactic background radiation field ionizes and heats the gas, we find that grain 
photoelectric heating by the FUV background is not sufficient to balance the C II 158 /im cooling rate 
inferred from damped Lya systems In these systems, a local energy source is required. We show that in 
the case of the z=1.919 damped Lya system toward Q2206— 19, the local source is FUV emission from 
the associated galaxy found by M0ller et al. (2002): the mean intensity inferred from photometry is in 
good agreement with the intensity J^ tars required to explain the cooling rate. The FUV mean intensity 
predicted for a cold neutral medium (CNM) model, J;f ars =(1.7t 2 'g) x 10" 18 ergs cm" 2 s" 1 Hz -1 sr" 1 
(95% c.l.), is the largest expected from our C II* study of 45 damped Lya systems. This may explain 
why this is the only confirmed damped Lya system yet detected in emission at z > 1.9. We argue that 
in most damped Lya systems with detected C II* absorption, J* tars is between 10~ 19 and 10~ 18 ergs 
cm" 2 s" 1 Hz -1 sr" 1 and heats the gas which is a CNM. By contrast, in most damped Lya systems 
with upper limits on C II* absorption the gas is a warm neutral medium (WNM). Surprisingly, the 
upper limits are compatible with the same range of Jf, tars values suggesting the majority of damped Lya 
systems are heated by radiation fields generated by a limited range of star formation rates per unit H I 
area, between 10~ 3 and 10 -2 . M Q yr" 1 kpc" 2 . We also show that C II* absorption is unlikely to arise 
in gas that is ionized. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 



1. INTRODUCTION 

External background radiation is the widely accepted 
heat source for the Lya forest. Defined to have H I col- 
umn densities 2V(H I) < 10 17 cm" 2 the Lya forest clouds 
are optically thin at the Lyman limit. As a result, ioniz- 
ing background radiation, i.e., radiation with photon en- 
ergies hv > hvn — 13.6 eV, penetrates the interior of 
the clouds, ionizing them such that the neutral gas frac- 
tion n(H°)/n<10 -5 , and heats them to temperatures in 
excess of 10000 K (Rauch 1998). Indeed, several authors 
have constructed scenarios within the context of CDM cos- 

iVisiting Astronomer, W.M. Keck Telescope. The Keck Obser- 
vatory is a joint facility of the University of California and the Cali- 
fornia Institute of Technology. 



mogonies in which most of the baryons in the universe at 
reside in the Lya forest and are photoionized by 
background radiation (e.g., Miralda-Escude et al. 1996; 
Machacek et al. 2000). 

Calculations by Haardt & Madau (1996; 2003) show 
that radiation emitted from the integrated population of 
QSOs and galaxies results in a mean background intensity 
J„ H ~ 10" 21 ergs cm" 2 s" 1 Hz" 1 sr" 1 at z ~ 3. Indepen- 
dent tests, such as the decrease with redshift of the number 
of absorbers per unit redshift in the vicinity of the QSO, 
i.e., "the proximity effect", indicate similar values for J VH 
(e.g., Scott et al. 2002). Ionizing background radiation 
has also been invoked to explain the ionization structure 
of metals in the Lya forest, and in higher column-density 
QSO absorption systems selected for C IV AA 1548, 1551 
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absorption (e.g. Mo & Miralda-Escude 1995), Lyman limit 
absorption (Prochaska 1999), and Mg II AA 2796, 2803 
absorption (e.g. Churchill, Vogt, & Charlton 2003). The 
predicted temperatures in these absorbers are also about 
10000 K and the neutral gas fraction, n(H°)/n << 1. 

In this paper we consider the nature of the heat source 
for another class of QSO absorption systems, the damped 
Lya systems (see Storrie-Lombardi & Wolfe 2000). While 
the damped Lya systems are subject to the same exter- 
nal radiation fields as the other absorbers, they differ in 
one important respect: the gas is neutral (e.g. Vladilo 
et al. 2001; Prochaska et al. 2002). Defined to have H 
I column densities N(H I)>2xl0 20 cm~ 2 , damped Lya 
systems have Lyman-limit optical depths, r VH >1000. As 
a result, the same photons that ionize and heat the Lya 
forest do not penetrate the interior of the damped Lya 
systems. Rather, these gas layers are penetrated only by 
background radiation with photon energies hv > 400 eV 
and hv < hvn- While the low-energy FUV photons do 
not ionize H, they do ionize neutral species of atoms with 
ionization potentials less than hvu such as C, Fe, and Zn, 
and more importantly, could contribute to the heat input 
of damped Lya systems through grain photoelectric heat- 
ing (Bakes & Tielens 1994). Similarly, X-ray photons can 
also heat the gas though photoionization. 

Does such background radiation dominate the heat in- 
put to the damped Lya systems? We shall address this 
question empirically since cooling rates have recently been 
measured for a representative sample of ~ 50 damped Lya 
systems (see Wolfe, Prochaska, & Gawiser 2003, hereafter 
referred to as WPG, for a discussion of the original sam- 
ple of 33 damped Lya systems). The answer has impor- 
tant implications. If heating by external sources does not 
match the measured cooling rates, then internal sources 
of heating must be found. This would distinguish damped 
Lya systems from all other classes of QSO absorption sys- 
tems, which are plausibly heated by background radiation. 

This paper is organized as follows. In § 2 we describe our 
technique for obtaining the thermal equilibria of damped 
Lya systems heated by external background radiation. In 
particular we compute the thermal equilibria of neutral gas 
layers exposed to background radiation fields recently com- 
puted by Haardt & Madau (2003). These authors com- 
puted backgrounds originating from the integrated popu- 
lation of (1) QSOs alone, and (2) QSOs and galaxies. To 
maximize the heating rates we consider those backgrounds 
that include galaxies because the galaxy contribution dom- 
inates at hv<hvH, where FUV radiation (6 < hv < 13.6 
cV) heats the gas by the grain photoelectric mechanism 
(WPG). This is the same process by which the neutral ISM 
of the Galaxy is heated (e.g. Bakes & Tielens 1994), and 
it is relevant here given the evidence for dust in damped 
Lya systems (Pettini et al. 1994; Pei & Fall 1995; WPG). 
At hv > 400 eV, soft X-rays heat the gas by photoionizing 
H, He, and abundant elements. X-ray heating occurs pri- 
marily through the photoionization of H and He because 
of the low metallicities of damped Lya systems. In §3 
we apply these results to test the background heating hy- 
pothesis. We compare the predicted heating rates to the 
158 /zm cooling rates deduced from the strength of C II* 
A 1335.7 and damped Lya absorption (see WPG). We 
study the multi-phase structure of the z = 1.919 damped 



Lya system toward Q2206— 19, and then discuss results 
for the full sample of 45 damped Lya systems. In § 4 
we examine the nature of possible internal heat sources 
by focusing on the 2=1.919 damped Lya system toward 
Q2206— 19. This damped Lya system is ideal for such a 
study because the presence of an associated galaxy at the 
same redshift, which is detected via its rest-frame FUV 
emission (M0ller et al. 2002), allows one to make an in- 
dependent estimate of the FUV radiation field incident on 
the absorbing gas. Since Fe is believed to be depleted in 
this DLA, the radiation is incident on dust, and the result 
is grain photoelectric heating. WPG show how to deduce 
the local FUV radiation field from the strength of C II* 
absorption. In § 5 we investigate the energy input for all 
the damped Lya systems in our sample, i.e., those with 
positive detections and those with upper limits on C II* 
absorption, to determine whether they are heated by back- 
ground radiation alone or require internal heat sources. In 
§ 6 we discuss the possibility that C II* absorption arises 
in ionized gas. Finally, concluding remarks are given in § 
7. 

Throughout this paper we adopt a cosmology consistent 
with the WMAP (Bennett et al. 2003) results, (Q. m , A , 
h) = (0.3, 0.7, 0.7). 

2. THERMAL EQUILIBRIA OF BACKGROUND HEATED 
DAMPED Lya SYSTEMS 

Consider a slab of neutral gas exposed to an isotropic 
background radiation field with mean intensity J v . We 
first compute the heating rates due to soft X-rays and 
FUV radiation. We then calculate the thermal equilibria 
resulting from a balance between heating and cooling. 

2.1. Heating 

Wolfire et al. (1995; hereafter W95) find the primary 
ionization rate of atomic species i due to soft X-rays to be 

Cxr = **j ^exph^A^HI)]^ s" 1 , (1) 

where <j l v is the photoionization cross-section per H atom 
of species i and the factor ct„V°(H I) is the optical depth 
due to H, He, and abundant metals such as C, O, Fe, etc. 
in an attenuating column density iV a (H I). To simulate 
conditions in the midplane of the gas layer we let 7V a (H 
I) = 0.5V(H I), where N(R I) is the H I column density 
detected along the line of sight. We compute the total ion- 
ization rate of H and He, including secondary ionizations, 
using the prescription of W95 including their modifica- 
tions of the results of Shull & van-Steenberg (1985). The 
resulting heating rate is given by 

Txr = 4vrY^ J ^-exp{-<j u N a (m)}<jSE h (E\x)dv , 

(2) 

where the summation is over species suffering primary ion- 
ization, and Eh(E l ,x) is the energy deposited as heat by 
a primary electron originating from atomic species i with 
energy E l in gas with electron fraction x. We used the re- 
sults of Shull & van-Steenberg (1985) to compute Eh, and 
an updated version of the Morrison & McCammon (1983) 
results to compute a l v and a v (McCammon 2003). 
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To compute the grain photoelectric heating rate due to 
FUV radiation we used the Weingartner & Draine (2001) 
expression 

r d = 10- 26 K Go/(G \/T/n e ,c 4 )exp(-T^) , (3) 

where Td is in units of ergs s _1 H^ 1 , k is the dust-to-gas 
ratio of the damped Lya system relative to that of the 
ISM (ISM refers to the ISM of the Milky Way Galaxy) , 
Go is 47T J J v dv integrated between 6 and 13.6 eV and is in 
units of 1.6xl0~ 3 ergs s _1 cm~ 2 , /, which is proportional 
to the heating efficiency, is a function of GoVT/n e and 
the fitting constants, Cj (where i=0— >6), and is the dust 
optical depth corresponding to N a (H. I). 

Weingartner & Draine (2001) show that Td is especially 
sensitive to two quantities: grain composition, i.e. the 
fraction of dust in silicate and carbonaceous grains, and 
the fraction of dust in small grains (radius < 15 A). Both 
are crucial because computations carried out so far show 
that small carbonaceous grains are the most efficient sources 
of photoelectric heating. The detection of the 2175 A ab- 
sorption feature in two damped Lya systems provide ev- 
idence for small carbonaceous grains at z < 1 (Junkkari- 
nen ct al. 2004; Motta et al. 2002), but the evidence 
is less clear for the redshift range of our sample. Pei, 
Fall, & Bechtold (1991) found no evidence for the 2175 
A feature, yet detection was expected based on their red- 
dening model. Thus damped Lya systems with z > 2 
could be comprised of silicate dust alone. However, all the 
published models for grain photoelectric heating include a 
population of small carbonaceous grains. 

To address this problem we consider models with a range 
of heating efficiencies that should include that of silicate 
dust. We adopt a model in which the fraction of car- 
bon locked in PAHs per H nucleus, 6c=0, and extinction 
ratio i?y=3.1, which is valid for the type of diffuse gas 
expected in damped Lya systems (see Table 2, row 19 
in Weingartner & Draine 2001); where Ry=Ay /Eb-v > 
Ay is the visual extinction, and Eb-v is the color ex- 
cess. We also assume a radiation spectrum approximat- 
ing Draine's (1978) fit to the interstellar radiation field in 
which case G « J„/(10 -19 ergs cm~ 2 s _1 Hz -1 sr _1 ) at A 
= 1500 A (see WPG). Application of this approximation to 
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WD- high 


[Si/H] 
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6xl(T 5 


3.1 
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[Si/H] 
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1.26xl(T 6 
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"Carbon abundance relative to solar. 
6 dust-to-gas ratio relative to the Milky Way ISM 
c Abundance of C atoms per interstellar H nucleus in PAHs 
d Ry= Ay /Eg_y, where Ay is visual extinction 

and Eb-v is the color excess. 
e (l) Weingartner & Draine 2001; (2) Bakes & Tielens 1994 
/ 10 [Si/H]( 1Q -0.2 _ 10 [Fo/Si]) 

91 [Si/H]( 1 _ 10 [Fc/Si]) 



background radiation involves no loss in generality since 
the background and interstellar radiation fields at hv < 
hvn arise from similar stellar populations (e.g. Bruzual & 
Chariot 2003). The resulting heating efficiency is the low- 
est predicted for the Weingartner & Draine (2001) models. 
The highest heating efficiency is predicted for their model 
with Ry=3A, &c=6xl0~ 5 , and an interstellar radiation 
spectrum (see their Table 2, last row). We shall examine 
the following three dust compositions: (1) a mix of sili- 
cate and carbonaceous grains in which &c= and i?y=3.1, 
(2) a mix of silicate and carbonaceous grains in which 
&c=6xl0~ 5 and Ry=3A, and (3) carbonaceous grains 
alone within the heating model of Bakes & Tielens (1994). 
In this way we shall test the sensitivity of our results 
to changes in grain properties. We show below that our 
results concerning the plausibility of background heating 
arc independent of the properties of the grains, and depend 
only on the carbon abundance and assumptions about 
thermal equilibrium. 

Next wc use the WPG method for computing n. Our 
prescription computes the fraction of Fe in grains and then 
assumes the number of depleted C or Si atoms per de- 
pleted Fe atom to be the same in damped Lya systems 
as in the Galaxy ISM. The depletion of Fc is calculated 
from the ratio of the observed amount of Fe to an un- 
depleted element: Zn or S would be ideal for this pur- 
pose. WPG use Si because (i) Si is measured over a wider 
range of rcdshifts than Zn, (ii) S is difficult to measure 
since the observed S II transitions occur frequently in the 
Lya forest, and (iii) though Si is undoubtedly depleted, 
the observed depletion level is very low; i.e., [Si/S]~0. 2 
This does not lead to inconsistencies in the case of sili- 
cate dust because the computed values of n are consis- 
tent with Si being depleted typically by less than 20 %. 

WPG find k = 10[ Si / H l"" (l0[ F °/ Si l-* - 10[ F °/ Si ls-) where 

[Si/H]j ni is the intrinsic abundance of Si relative to H, 
[Fe/Si]j n £ is the intrinsic abundance of Fe relative to Si, 
and [Fe/Si] gas is the measured abundance of Fe relative 
to Si. Following WPG we assume minimal depletion oc- 
curs when [Fe/Si]i n t=— 0.2, i.e., an intrinsic enhancement 
of a elements such as Si, and maximum depletion occurs 
when [Fe/Si]j n t=0. Therefore, we adopt two Weingartner 
& Draine (2001) models, which should bracket a plausible 
range in heating rates. The "WD-low" model combines 
dust composition (1) with minimal depletion: this model 
is similar to one of the "SMC" models considered by WPG 
and results in the lowest heating rates. 3 The "WD-high" 
model combines dust composition (2) with maximal de- 
pletion: this model produces the highest heating rates. 
For completeness we combine the Bakes & Tielens (1994) 
model of carbonaceous grains with the assumption of max- 
imal depletion: this model, referred to as the "BT" model, 
produces heating rates between the predictions of the other 
two models, and is one of the "Gal" models considered by 
WPG. The properties of the models are summarized in 
Table 1. 

The background mean intensities at redshift z, J v {z), 

2 The abundance ratio of elements X and Y is given by 
[X/Y] =logi (X/Y)-logi (X/Y) s . 

^Minimal depletion does not apply to damped Lya systems in 
which [Fe/S i]gas > — 0.2. For these objects the "WD low" model is 
a hybrid of 6c =0 heating rates and maximal depletion. 



4 



STAR FORMATION IN DAMPED Lya SYSTEMS 



were computed with the software package CUBA written 
and made available to us by Haardt & Madau (2003). 
These authors give a detailed description of how J v (z) 
is calculated from the luminosity function and spectra of 
QSOs and how the radiation is reprocessed by the inter- 
vening gas (Haardt & Madau 1996). Because grain pho- 
toelectric heating depends on J v at photon energies, hv < 
hvH, where galaxies dominate QSOs, the galaxy contribu- 
tion to the background radiation is crucial. For this rea- 
son we considered backgrounds that include the integrated 
contribution of galaxies. We used models in which the ra- 
diation escaping from the galaxies is attenuated by dust 
because this is the radiation that contributes to the back- 
ground intensity incident on a given damped Lya system. 
In this case the inferred luminosity per unit bandwidth 
per unit comoving volume decreases with redshift at z > 
1.5, as found by Giavalisco et al. (2003). The predicted 
background spectra at the redshifts of five damped Lya 
systems from our sample are shown in Figure 1. The fac- 
tor of ~ 60 increase in J v as v decreases across vh is the 
most striking feature in Figure 1. This feature is missing 
from backgrounds computed from QSOs alone (see Haardt 
& Madau 1996) and arises from stellar radiation, which is 
optically thick at the Lyman limit, and from an assumed 
escape fraction from galaxies of lO^ 1 at hv > hvn (Haardt 
& Madau 2003). Therefore, J v at hv < Kvr is typically 
10~ 20 ergs cm~ 2 s _1 Hz -1 sr _1 , which is a factor of ten 
higher than J v produced by QSOs alone. By comparison, 
the ambient FUV radiation in the ISM is characterized by 
J v <~ 10~ 19 ergs cm~ 2 s _1 Hz -1 sr _1 . The spectral regions 
that are critical for DLA heating correspond to 6 < hv < 
13.6 eV and hv > 400 cV. Figure 1 shows that for a given 
dust-to-gas ratio and heating efficiency the FUV heating 
rates should decrease no more than 0.5 dex between z=1.6 
and 4.3, while for a given N(H I) and metallicity the X-ray 
heating rates should increase by no more than 0.3 dex in 
the same range of redshifts. 

2.2. Cooling 

Cooling rates were computed according to the prescrip- 
tion given by W95 and WPG. Thus, we include cooling 
due to Lya emission, and grain radiative recombination, 
which dominate at high temperatures, T > 3000 K. At 
lower temperatures, excitation of fine-structure states of 
abundant ions such as C + , O , Si + , and Fc+, and excita- 
tion of metastable states of O dominate the C II cooling 
rate. In all cases, [C II] 158 /im emission, which arises 
through the transition between the 2 Pz/2 an d 2 -f\/2 fine- 
structure states in the 2s 2 2p term of C + , is the domi- 
nant coolant in the CNM (cold-neutral medium) phase, 
while Lya emission dominates in the WNM (warm-neutral 
medium) phase. Because CMB excitation of the 2 Ps/2 
state in C II increases with redshift, the spontaneous emis- 
sion rate of 158 /um radiation per H atom, l c , need not 
equal the cooling rate in low density gas at high redshifts. 
To separate the effects of cooling from radiative excita- 
tion we note that l c =nAcn+(l c )cMB, where nAcn is the 
cooling rate due to [C II] 158 emission: the cooling rate is 
defined as the difference between collisional excitation and 
de-excitation rates of the 2 Pz/2 state times the excitation 
energy, hv u i, of the 2 P^/2^ 2 P\/2 transition. The quantity 
{I c )cmb is the spontaneous emission rate of 
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Fig. 1. — Haardt & Madau (2003) background spectra for damped 
Lya systems with redshifts spanning the redshift interval of our 
sample. Spectral features are explained in the text 



158 /im radiation per H atom in the limit of pure radia- 
tive excitation by the CMB, and equals 2(C/H)A u ihv u i 



exp 



-hv u i/[k(l + z)T C mb] 



where (C/H) is the carbon 



abundance. The quantity A u i is the emission coefficient 
for spontaneous photon decay, and the current tempera- 
ture of the CMB is T CM b=2.728 K. 

3. HEATING BY EXTERNAL BACKGROUND RADIATION 

In this section we apply the formalism of § 2 to study 
the thermal structure of damped Lya systems subject to 
external background radiation. 

3.1. Background Heating of the z = 1.919 Damped Lya 
System Toward Q2206— 19 

We tested the background heating hypothesis by com- 
paring the predicted values of l c with observations. The 
observed rate, l° bs = N(C ll*)A ul hv ul /N(R I) where N(C 
II*) is the observed column density of C II in the excited 
2 P 3 / 2 fine-structure state. The C II* column densities were 
derived from accurate measurements of C II* A 1335.7 ab- 
sorption lines arising from the 2 P^/2 state. WPG deter- 
mined l° bs for 33 damped Lya systems. Most of the data 
were acquired with the HIRES echelle spectrograph (Vogt 
et al. 1994) on the Keck I 10 m telescope, and in three ob- 
jects with the UVES echelle spectrograph on the VLT 8 m 
telescope. More recently, l° bs was obtained by Prochaska 
et al. (2003) from measurements for 16 damped Lya sys- 
tems with the ESI echelle spectrograph imager (Sheinis et 
al. 2002) and with UVES for another system (Pettini et 
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Fig. 2. — Two-phase diagrams for gas in the 2=1.919 damped Lya system toward Q2206— 19. The heat input is FUV and soft X- 
ray background radiation. The heating mechanisms are grain photoelectric emission for FUV radiation and photoionization of H, He, and 
abundant elements for the X-rays. In this case we assume the "WD-low" model and minimal depletion (see text). The relevant data for this 
damped Lya system are in Table 2. Panel (a) shows pressure versus density. Panel (b) shows fractional ionization, x, and electron density, 
n e versus density. Magenta curves in panel (c) show grain photoelectric heating (PE), X-ray heating (XR), and C I photoionization heating 
rates versus density. Dotted blue curve is [C II] 158 jim fine-structure cooling rate, and green curves are Lya and grain recombination cooling 
rates. Dotted red curve is [S II] 35 fim fine-structure cooling rate. Solid blue curve is [C II] 158 fim spontaneous energy emission rate. Black 
curve (r) is total heating rate. Black horizontal solid line is logioi" 68 and black dotted horizontal lines are 2-cr errors in logioZ" 6 " Panel (d) 
shows temperature and C II/C I ratio versus density. 



al. 2002). Thus, the total sample consists of l° bs measure- 
ments for some 50 damped Lya systems. In this paper we 
shall analyze the properties of 45 of these objects. These 
include 23 damped Lya systems with positive detections 
of C II* absorption and 22 with upper limits on C II* 
absorption. 

To predict l c we computed thermal equilibria resulting 
from the energy balance condition, T—nA. The total heat- 
ing rate T=TxR+^d, where the radiation input is from 
XR and FUV background radiation alone; i.e., no inter- 
nal heat sources are included. The quantities n and A are 
gas volume density of hydrogen and total cooling function 
respectively. Thermal equilibrium constrains the tempera- 
ture, pressure, cooling rates, etc. to be unique functions of 
density. It also results in l c as a unique function of density, 
i.e., l c (n). 

In order to interpret the results for our full sample of 45 
damped Lya systems the reader is referred to Figure 2 for 
the phase diagrams of the z = 1.919 damped Lya system 
toward Q2206— 19, one of our sample damped Lya sys- 
tems, which is hereafter referred to as the Q2206A damped 



Lya system. The figure shows equilibrium curves for (a) 
pressure, (b) electron fraction and electron density, (c) 
cooling rates, heating rates, and l c , and (d) temperature 
and C II/C I column-density ratios (hereafter C II/C I de- 
notes the ratio N(C U)/N(C I)) as functions of hydrogen 
density. The equilibria are computed for the "WD-low" 
model. The curves show the classical two-phase equilib- 
ria (see Field, Goldsmith, & Habing 1969; W95) exhibited 
by neutral gas layers heated by the grain photoelectric 
mechanism and ionized by soft X-rays (W95, WPG). In 
that case the equilibria admit the presence of a thermally 
stable WNM at low densities and high temperatures sepa- 
rated from a thermally stable CNM at high densities and 
low temperatures by a thermally unstable intermediate- 
density region in which d(logP) / d(logn) < 0. 4 Figure 2c 
shows nAcn (dotted blue curve) to increase with density. 

4 The thermally unstable region occurs at lower densities than in 
the ISM of the Galaxy because the FUV heating rate produced by 
external background radiation incident on gas with a low dust-to-gas 
ratio is about 100 times lower than in the ISM where the heat source 
is internal FUV radiation emitted by stars in the Galaxy (see WPG). 
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The increase is steep at n < 10" 2 cm" 3 , very steep at n w 
10~ 2 cm" 3 , flatter at 10 -2 < n < 10" 1 cm" 3 , insensitive 
to n at 10" 1 < n < 10 +2 cm" 3 , and steeper again at n > 
10+ 2 cm" 3 . 

To understand this behavior we note that the cooling 
rate is dominated by processes other than 158 /zm emission 
at n < 10" 1 cm" 3 . Consequently, nAcn does not track the 
total heating rate, T, which is flat at the latter densities. 
The explicit form of nAcn under these conditions is given 
by 

nAcn cx n e (C/H)r!( 2 P 3/2 , 2 J P 1/2 )T- 1 /2 ex p(_^/A : T) , 

(4) 

where the collision strength, f2( 2 P 3 / 2 , 2 ^1/2)1 is relatively 
insensitive to changes in n and T (see Blum & Pradhan 
1992). At n < 10" 2 cm" 3 the increase in nAcn is driven 
primarily by the increase in electron density with increas- 
ing n (see Figure 2b) because (1) the Boltzmann factor, 
exp(—hv u i/kT), is constant since T » hv u ijk~ 92 K, 
and (2) T is insensitive to n at these densities (see Fig- 
ure 2d). At n w 10" 2 cm" 3 the very steep drop in T 
marks the transition between WNM and CNM gas in the 
region of thermal instability due to the shift from Lya 
cooling to fine-structure line cooling by low ions of abun- 
dant elements. The sharp increase in nAcn is due to the 
increase in the T" 1 / 2 factor because at these tempera- 
tures the Boltzmann factor is constant. The decrease in 
T also increases the radiative recombination coefficients, 
which cause the abrupt drops seen in x, n e , and C II/ C 
I. When 10" 2 < n < 10" 1 cm" 3 the increase in nAcn 
is modulated by the Boltzmann factor which decreases as 
T becomes comparable to hv u i/k. However, at n > 10" 1 
cm" 3 the dependence of nAcn on n is dictated by the 
density dependence of the total heating rate, r(n), since 
nAcn dominates the cooling rate at these densities and 
thus equals the total heating rate. This accounts for the 
flat behavior of nAcn at 10" 1 < n < 10 +2 cm" 3 where 
the grain photoelectric heating rate dominates and is in- 
dependent of density, and the steep increase of nAcn at 
n > 10 2 cm" 3 where the C I photoionization heating rate 
per atom becomes dominant. Note that l c (n) departs from 
nAcn only at n < 10" 2 cm" 3 where radiative excitations 
by CMB radiation dominate collisional excitations. 

Comparison between predicted and observed values of 
l c in the case of the Q2206A damped Lya system shows 
that heating and ionization by background radiation can- 
not account for the observed 158 ^tm cooling rate unless 
the densities are higher than n > 10 3,4 cm" 3 for the "WD- 
low" model shown in Figure 2 (specifically see Figure 2c). 
But in the case of FUV background heating alone , densi- 
ties this large are ruled out by the observational constraint, 
C II/C I > 2x 10 4 . The connection between n and C II/C I 
follows from the condition of photoionization equilibrium, 
which implies 

en Go 

C I K a{C\)xn ' [ ' 

where a(CI) is the C I radiative recombination coefficient. 
In thermal equilibrium C II/C I and x are functions of 
n, which are plotted in Figure 2. The low value of Go 
(=0.14) and the large lower limit on C II/C I combine to 
set an upper limit on density, n < 10" 2 cm" 3 for FUV 
background heating alone. From the l c (n) curve in Figure 



2c we see this corresponds to l c < 10~ 278 ergs s" 1 H" 1 
which should be compared to l° bs > 10" 26 5 ergs s" 1 H" 1 ; 
i.e., the 2— a lower limit denoted by the lower dotted hor- 
izontal line. This demonstrates that background heating 
alone cannot account for the value of l° bs . 

Although this discrepancy was derived for the "WD- 
low" grain model, it is quite general as it is insensitive 
to grain composition. To see this we compare l c {n) with 
l° bs for all 45 damped Lya systems in Table 2. In each 
damped Lya system we compare the predictions of the 
"WD-low" model with the "WD-high" and "BT" mod- 
els. Before discussing the full sample we focus first on the 
Q2206A damped Lya system. Figure 3a (first column 
and second row from the top) shows the l c {n) predicted 
by each model is more than a factor of 10 below the lower 
limit on l° bs at densities, n<10" 2 cm" 3 , permitted by the 
C II/C I constraint. The discrepancies between l c (n) and 
l obs differ slightly for each model for reasons that are inde- 
pendent of grain properties. First, notice that at n <10" 2 
cm" 3 the Z c (n) predicted by the "WD-high" and "BT" 
models are degenerate and higher than predicted by the 
"WD-low" model. At such low densities the gas is a WNM 
where Z c (n) oc n e (C/H). Because there are no unsaturated 
transitions of C II we assume C traces Fe and as a re- 
sult we let [C/H]=[Si/H]+[Fe/Si]j„ t to determine the car- 
bon abundance (see § 5.1 in WPG). From Table 1 we see 
that the carbon abundance, [C/H], is lower in the "WD- 
low" model than in the "WD-high" or "BT" models be- 
cause [Fe/Si]j n t=— 0.2 for the minimal depletion assumed 
in the "WD-low" model and [Fe/Si]j„t=0 for the maximal 
depletion assumed in the "WD-high" and "BT" models. 
Therefore, l c (n) is predicted to be lower in the "WD-low" 
model. Second, the C II/C I constraints on l c (n) will differ 
because the observational limits on C II/C I will be lower 
by 0.2 dex in the "WD-low" model. This is because N{C 
II) will be 0.2 dex lower while N(C I) will be the same, 
as it is inferred directly from the strength of C I A 1656.9 
absorption. The predicted dependence of C II/C I on n 
will also differ in each model because the T(n) and n e (n) 
relations differ. The important point is that the C II/C I 
limits restrict gas heated by background radiation to the 
WNM phase where the limits on l c (n) depend only on the 
carbon abundance and the assumption of thermal equilib- 
rium, both of which are independent of grain properties. 
This is a generic trait of the full sample to which we now 
turn. 

3.2. Background Heating of the Full Sample 

We repeated these procedures for the 44 remaining damped 
Lya systems for which it was possible to calculate thermal 
equilibria for all three dust models. This yielded phase di- 
agrams for a total of 23 damped Lya systems in which C 
II* A 1335.7 was detected and 22 with upper limits. Com- 
parisons between l° bs and the predicted l c {n) curves are 
plotted for the detections in Figures 3a and 3b, and for 
the upper limits in Figures 4a and 4b. Relevant data for 
all of the damped Lya systems are given in Table 2. The 
coordinate name of the background QSO is given in col- 
umn 1, the absorption rcdshift in column 2, the value of 
logmC II/C I in column 3, the observed [Fe/H] and [Si/H] 
abundance ratios in columns 4 and 5, \ogu)l° bs in column 6, 
and the mean intensity due to starlight, Jf ars 7 in column 
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Fig. 3. — l c {n) curves for 12 damped Lya systems with positive 
detections of C II* A1335.7 absorption. The l c (n) are computed 
from thermal equilibria where heat input is from background radia- 
tion alone. Solid curves show l c (n) for the "WD-low" model, dashed 
curves for the "BT" model, and dotted curves for the "WD-high" 
model. In all cases the heating mechanisms are grain photoelectric 
heating and X-ray photoionization. Data are shown as solid hori- 
zontal lines for logio?° 6s and dotted lines for 1— a errors in \ogu)l° bs . 
The solid horizontal lines in Q2206-19 and Q0458-02 damped Lya 
systems represent 2— a lower limits on logio(° i>s . Arrows on l c (n) 
curves are upper limits on n imposed by observed lower limit on 
C II/C I. These become data points and error bars for Q1232+08, 
Q1331+17, Q2231-00, and Q0812+32 where C II/C I is measured. 



7 (see § 5). 

In Figures 3a and 3b we compare the predicted l c (n) 
with l° bs for damped Lya systems with positive detec- 
tions. In every case in Figure 3a and in 9 out of 12 cases 
in Figure 3b the predictions fall significantly (more than 
4c; ogi0 / e ) below the observations at densities permitted by 
the C II/C I ratios (shown as horizontal arrows directed 
toward decreasing n). Most of the predicted l c (n) curves 
exhibit behavior similar to that described for the Q2206A 
damped Lya system. Specifically, the "WD-high" and 
"BT" models are degenerate at low n, but diverge at n 
sufficiently large that l c (n) equals the grain photoelectric 
heating rate. At the latter densities l c {n) is higher for 
the "WD-high" model than the "BT" model because of 
the difference in heating efficiencies. Furthermore, l c {n) 
for both "WD-high" and "BT" models are higher than 
for the "WD-low" models because the heating efficiencies 
and k (Table 1) are higher in the "WD-high" and "BT" 
models. However, there are differences with the Q2206A 
result. First Figure 3b shows that in some damped Lya 
systems with very low dust-to-gas ratios (e.g. Q0953+47 
and Q0747+27) l c (n) decreases rather than increases with 
n above densities corresponding to the WNM— CNM tran 
sition zone. This behavior stems from the low values of k, 
which suppress grain photo-electric heating in favor of 
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X-ray heating, which decreases with n (see Figure 2c). 
Second, l c (n) exhibits a flat plateau at low densities in 
damped Lya systems with large redshifts since l c (n) ap- 
proaches (I c )cmb in the low density limit. Because {I c )cmb 
« l° c bs , the CMB has only marginal effect on the [C II] 
158 (im emission rate for every damped Lya system in 
Figure 3a and most damped Lya systems in Figure 3b. 

On the other hand, background radiation is the likely 
source of (/ c )o6s f° r two of the 11 positive detections in 
Figure 3b. Specifically, the CMB is the dominant factor 
for the Q0255+00 and Q2241+13 damped Lya systems 
where (I c )cmb ~ {l c )obs- It accounts for l° bs at all den- 
sities in all three grain models for Q2241+13 5 , and it 
accounts for l° bs in the "WD-high" and "BT" models for 
Q0255+00. But since the "WD-low" model predicts l c (n) 
to be 2.6ai ogio i c below l° bs at low n, a heat source in ad- 
dition to background radiation may be required for the 
latter damped Lya system. The CMB is unlikely to ex- 
plain l° bs for any of the remaining damped Lya systems 
with positive detections of C II* absorption. Therefore, 
with one exception we find that \ogn)l° bs — logio/ c (") > 
3czogi ic & t densities permitted by the C II/C I constraint 
when l c {n) is calculated under the assumption that exter- 
nal background radiation is the sole source of heating. The 
exception is the Q 1346— 03 damped Lya system (see Fig- 
ure 3b) in which heating by X-ray background radiation 
can account for l° bs when 1 < n < 25 cm -3 . Therefore, we 
conclude that background radiation alone cannot explain 

5 The "WD-low" model for the Q2241+13 damped Lya system is 
degenerate with the "WD-high" and "BT" models in Figure 3b since 
we assumed maximal rather than minimal depletion. This occurs 
because [Fe/Si] 9as > - 0.2 in this and in the Q0201+11, Q1021+30, 
Q1506+52, and Q2344+12 damped Lya systems. 
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"Data from WPG 

'Data presented here for the first time 

c Estimates derived assuming [C/H] = [Si/H], which is appropriate for the "WD-high" and "BT" model. Estimates for the "WD- low" models arc 
0.2 dex lower (see Table 1 

d In this case errors are 2— a 

e Wavclcngth regions covering C I transitions not observed 

•^Entries are CNM solutions for J^ tars ( in units of 10 -19 ergs cm -2 s _1 Hz -1 sr _1 ) in damped Lya systems with positive detections. Numbers 
in parentheses denote WNM solutions for cases where WNM hypothesis is plausible. 



the inferred 158 emission in 20 out of the 23 positive 
detections in Figures 3a and 3b. In each of the 20 cases 
an internal source of heating is required. 

Figures 4a and 4b show l c (n) curves predicted for the 22 
damped Lya systems with upper limits on l° bs . As for the 
damped Lya systems with positive detections, l° bs exceeds 
l c (n) in most cases. But there is an important difference: 
since l° bs in Figures 4a and 4b are upper limits, the dis- 
crepancy with the model predictions need not require the 
presence of internal heat sources. Thus, if the true l° hs re- 
vealed by future detections of C II* absorption intersected 
the l c {n) curves, then background heating could account 
for the C II cooling rates, provided the gas density sat- 
isfied the condition l c (n)=l° bs , which defines the density 
n=riBKD- Because ubkd lies in the WNM phase for every 
damped Lya system shown in Figure 4, the gas would be 
a WNM. Most of the damped Lya systems illustrated in 
Figure 4 are candidates for this condition since the upper 
limits on l° bs are typically less than 0.2 dex above l c (n). 
But, internal heating would be required if the true l° bs 



remain above the l c (n) curves for all n. This is clearly 
possible for the damped Lya systems toward Q1223+17, 
Q2348-14, Q2206-19A, and Q0209+52 where the upper 
limits on l° bs are typically more than 0.5 dex above the 
peaks in the l c (n) curves. In § 5.2 we show that internal 
heating may also be present even when l° bs dips below the 
peaks of the l c (n) curves, which we emphasize are pre- 
dicted for the case of background heating alone. 

To summarize, the evidence discussed so far suggests the 
following: Most damped Lya systems with positive det 
ections of C II* absorption require internal heat sources, 
and we are able to rule out the hypothesis that damped 
Lya systems with upper limits on C II* absorption consist 
of CNM gas heated by background radiation alone. In 
§ 4 we discuss independent evidence for internal heating 
by starlight for one damped Lya system in which C II* 
absorption is detected. In § 5 we consider this possibility 
for all the damped Lya systems in Table 2, including 
those with upper limits on C II* absorption. 
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Fig. 4. — l c {n) curves for 12 damped Lya systems with upper 
limits on C II* A 1335.7 absorption. The i c (n) curves are computed 
as in Figures 3a and 3b. Resulting 2— a upper limits on logioi° i>s 
are solid horizontal lines with downward pointing arrows. 
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4. INDEPENDENT EVIDENCE FOR HEATING BY 
STARLIGHT IN A DAMPED Lya SYSTEM 

In § 3.2 we concluded that heating by external back- 
ground radiation is insufficient to account for the [C II] 158 
(im emission rates inferred for most damped Lya systems 
with detected C II* A 1335.7 absorption. Rather, an in- 
ternal heat source is required. WPG and WGP suggested 
grain photoelectric heating by FUV radiation from stars 
associated with the absorbing gas as the heat source. They 
neglected background radiation because the computed ra- 
diation fields available to them included the contribution 
by QSOs alone (e.g. Haardt & Madau 1996). These fields 
have a factor of 10 lower intensity than the galaxy contri- 
bution at hv < hvu and, therefore, make an insignificant 
contribution to the required heating rate. 

Evidence supporting grain photoelectric heating was pre- 
sented in WGP. In particular the statistically significant 
correlation between l° bs and dust-to-gas ratio, k, argues 
against mechanical heating mechanisms such as turbu- 
lent dissipation (e.g. Wolfire et al. 2003) and is instead 
strong evidence for grain photoelectric heating, provided 
l obs represents the heating rate. Accurate upper limits 
on Si II* A 1264.7 absorption established by Howk et al. 
(2004a, b) for three damped Lya systems with detected C 
II* A 1335.7 absorption give independent evidence that 
the C II* absorption arises in CNM gas and, as a re- 
sult, l° bs must equal the heating rate (see § 3.1). There- 
fore, the evidence accumulated so far indicates (1) pho- 
toejection of electrons from grains as the heating mech- 
anism, and (2) that background radiation does not pro- 
vide sufficient photons. While this points toward locally 
produced FUV radiation, i.e. starlight, as the source of 
photons, WGP did not offer direct evidence for a link be- 
tween the inferred radiation intensity and starlight in a 
given damped Lya system. 

In this section we present evidence for such a link. Specif- 
ically, M0ller et al. (2002) detected starlight emitted by a 
galaxy associated with the Q2206A damped Lya system. 
Therefore, it is possible to deduce J v directly and compare 
it with the J v inferred from 

C II* absorption at hv = 8 eV. That is the purpose of this 
section. 

4.1. Jf ars Inferred From C II* Absorption 

Figure 5a shows the C II* A 1335.7 absorption profile 
for the Q2206A damped Lya system. This spectrum was 
acquired with UVES, the echelle spectrograph on the VLT 
8 m telescope, by Pettini et al. (2002) and was obtained 
by us from the UVES archive. This transition is clearly 
detected between v=— 20 and v=+30 km s _1 . At v >+ 
30 km s _1 the C II* transition is blended with a broad 
Lya forest feature extending between v=+30 and v=+300 
km s _1 . If one assumes that C II* A 1335.7 cuts off at 
v=+ 30 km s^ 1 , then logi iV(C IP)=13.65±0.05 cm" 2 . 
This is the column density used to compute l° bs in Figure 
3. However, this is likely to be a lower limit to the actual 
column density. The reasons are illustrated in Figures 5b 
and 5c, which show the Si II A 1808.0 and Ni II A 1370.1 
velocity profiles (Prochaska & Wolfe 1997a). These and all 
other low-ion resonance-line profiles in this damped Lya 
system exhibit significant absorption between v=— 20 and 
u=+120 km s" 1 . In fact Prochaska & Wolfe (1997b) used 
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Fig. 5. — Panel (a) shows velocity profile of C II* A 1335.7 ab- 
sorption (solid curve) in the Q2206A damped Lya system. Dotted 
curve at lower velocities is C II A 1334.5 absorption. Dotted curve at 
higher velocities is a Lya forest line. Data acquired with UVES on 
the VLT. Panels (b) and (c) shows Si II A 1808.0 and Ni II A 1370.1 
velocity profiles in the same absorption system. Data acquired by 
Prochaska & Wolfe (1997a). Vertical line shows coincidence of the 
lowest velocity component in all three transitions. 



the low-ion profiles to obtain an absorption velocity inter- 
val, Au=130 km s _1 for the Q2206A damped Lya sys- 
tem. The C II* profile is likely to exhibit a similar velocity 
structure since the C II* and low-ion resonance-line veloc- 
ity profiles are indistinguishable in all but one damped 
Lya system for which C II* A 1335.7 is detected (sec 
WGP). When we integrate the Si II and Ni II profiles over 
the entire absorption range, we find the column densities 
to be 0.50 dex higher than when integrated from v=— 20 
to +30 km s _1 . Correcting the C II* column density 
by this factor we find logi iV(C H*)=14.15±0.10 cm" 2 . 
We take the latter value as an upper limit on N(C II*) 
and logm N(C II*) = 13.65 cm" 2 for the lower limit. By 
adopting the mean of these for the detected value, we find 
logi iV(C II*)=13.97^;^ where the errors are the differ- 
ences between the limits and the mean. The corresponding 
158 /im emission rates are given by logio^° fcs =— 26.20^q 32 
ergs s" 1 H _1 . We believe the errors in logio-/V(C II*) and 
\ogiol° bs correspond to 95% confidence intervals. 

To obtain the FUV mean intensity generated by lo- 
cal stars from l° bs we recompute the thermal equilibrium 
of the neutral gas. First we compute the grain photo- 
electric heating rate, Td, due to the total mean inten- 
sity, J„=Jf ars +J bkd , where Jf ars and J bkd are the FUV 
intensities due to local stars and background radiation. 
Then we recompute the X-ray heating rate by adding the 
contribution from hot gas assumed to be associated with 
local stars to the X-ray background intensity (see Wolfire 
et al. 1995 and WPG). We also include heating by cos- 



mic rays generated by local stars. In both cases we let 
the heating rates be proportional to the source function 
for J° tars ; i.e., the SFR per unit area, ip* (see WPG). 
We compute the total heating rates for a grid of mean 



intensities in the range 5x10 



-21 



< J' 



< 5x10 



-17 



ergs cm~ 2 s _i Hz -1 sr . We use the prescription dis- 
cussed in § 2 to obtain metal abundances and dust-to- 
gas ratios, which are used to generate cooling rates and 
heating rates respectively. Thermal equilibrium models 
give P(n) and l c (n) curves similar to those shown in Fig- 
ure 2. To obtain Jf, tars corresponding to a given l° bs we 
assume damped Lya systems to be two-phase media in 
which CNM gas is in pressure equilibrium with WNM gas 
at pressure Peq = {PminPmax)^ > where P m in and P m ax 
are the pressure minimum and maximum exhibited by a 
two-phase medium in thermal equilibrium (see Figure 2a). 
This results in two unique solutions for each l° bs , one in 
which C II* absorption arises in the CNM and the other 
in the WNM. In WGP we argued against the WNM hy- 
pothesis as it gives rise to more background radiation than 
observed at the current epoch. The argument against the 
WNM hypothesis is further supported by the absence of 
Si II* A 1264.7 absorption in three damped Lya systems 
(see Howk et al. 2004a, b). Nevertheless, we consider the 
WNM solutions anew since they provide independent tests 
of the WNM hypothesis. 

Figure 6 shows the Q2206A phase diagrams correspond- 
ing to the CNM solution. In this case the mean inten- 
sity resulting in l c (n)=l° bs at the CNM density, iicnm, 
is given by J** ars =2.3x 10~ 18 ergs cm~ 2 s _1 Hz -1 sr -1 , 
where ticnm is the higher density root of the polynomial, 
P(n) — P eq =0. We find ncNM=l-3 cm~ 3 as shown by the 
vertical dot-dashed line which crosses the intersections be- 
tween l c (n) with l° bs and P(n) with P eq (Figures 6a and 
6c). We again adopt the "WD-low" model. 

Comparison between Figures 2 and 6 reveal several dif- 
ferences. First, P m ax and P m i n occur at higher densities 
when local radiation fields with intensities much higher 
than the background are included. These pressures oc- 
cur where fine-structure cooling dominates resonance line 
cooling and signify the transition from WNM to CNM gas. 
In the CNM gas the temperature is determined by the 
balance between [C II] 158 fim. cooling and grain photo- 
electric heating in the low temperature limit. Because the 
density threshold above which grain photoelectric heating 
dominates X-ray and cosmic ray heating is proportional to 
jstars^ transition to CNM occurs at higher densities as 
jstars increases (compare Figures 2c and 6c). The larger 
value of J„ also results in larger C II/C I ratios at a given 
n (compare Figures 2d and 6d). Therefore, the high den- 
sity solutions ruled out for background heating alone are 
allowed when local sources are included. Indeed the value 
of ficNM is consistent with the lower limit set on C II/C I 
for this damped Lya system. 

To make realistic predictions of Jf, tars for the "WD- 
low" model we proceeded as follows. First, we computed 
jstars corresponding to l° bs . We repeated the calcula- 

6 In WPG we generated the J„ grid starting with a grid of logioV>* 
between —4.0 and 0.0 Mq yr _1 kpc~ 2 . In this paper, we work with 
J„ directly since it can be inferred from observations. We consider 
the SFRs later in § 6 when we connect J£ tars to the rate of star 
formation. 
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Fig. 6. — Two-phase diagrams for the z = 1.919 damped Lya system toward Q2206— 19. Same as Figure 2 except internal sources of 
heating are included. Specifically, grain photoelectric heating by starlight with FUV mean intensity, J^ tars =1.7x 10 — 18 ergs cm -2 s — 1 Hz -1 
sr _1 is included. Heat inputs from locally generated soft X-rays and cosmic rays are also included (see text). This CNM solution yields 
agreement between l c (n) and l° bs at the CNM density, ncNM=l-3 cm -3 . This is the thermally stable density in which CNM and WNM gas 
are in pressure equilibrium at pressure, P=(PminPmax y/ 2 (see text). 



tion for the 95 % confidence lower and upper limits on 
l° bs . We then assumed the errors in Jl tars to be the dif- 
ferences between the limiting values of J s J ars and that 
corresponding to the mean l obs . These errors are due 
to systematic uncertainties in l° bs , and are larger than 
the random statistical errors for this damped Lya sys- 
tem. Other systematic errors arise from the uncertainty 
in P eq , k, and J bkd . To evaluate these we recomputed 
jstars by ^ i ettmg p eq vary between P min and P max , (2) 
varying k between the conditions of minimal and maxi- 
mal depletion, and (3) varying J bkd between 0.1 and 3.0 
times the Haardt-Madau backgrounds used here, which 
are overestimates of the uncertainties in J bkd (Haardt & 
Madau 2003). Again assuming the variances are deter- 
mined by the differences between these limiting values and 
the mean Jf ars defined by P eq = (P m i n Pmax) 1/2 , we find 
that J^ tars = (1.7+^)xl0- 18 ergs cm" 2 s" 1 Hz" 1 sr" 1 for 
the CNM solution and J^ ars =(3.9l^)x 10" 17 ergs cm" 2 
s _1 Hz -1 sr" 1 for the WNM solution. The latter two are 
the means of the minimal and maximal depletion "WD- 
low" solutions and the errors correspond to 95 % confi- 
dence intervals. In every case the J* tars predicted by the 
WNM solutions exceed the 95 % confidence upper limit 



for our CNM solution. These estimates for the "WD- 
low" model are also consistent with the predictions of most 
"WD-high" and "BT" solutions. Figure 7 plots the CNM 
solutions of Jf ars versus J bkd / \J bkd ) h m , where {J bkd ) H M 
is the Haardt-Madau background solution adopted here. 
We computed solutions for all the three grain models used 
in this paper. Comparison with the above estimate of 
jstars £ or fog CNM model shows consistency between the 
"WD-high", "BT", and "WD-low" models with the possi- 
ble exception of the "BT" model, which predicts J^, tars to 
be slightly lower than our estimate. 7 

Finally, it is interesting that these are the highest in- 

7 Note, J^ tars predicted by the "WD-high" model is larger than 
predicted by the "WD-low" model. This seemingly contradicts the 
higher heating efficiency of the "WD-high" model, because Jjf ars 
required to produce a given heating rate decreases with increasing 
heating efficiency. However, that argument holds for a given density. 
Because of the high metallicity of the Q2206A damped Lya system 
and the shape of the heating function, r,j(n), the value of iicnm 
is lower for the "WD-high" than for the "WD-low" model. As the 
density decreases, J° tars increases in order to balance the cooling 
rate, l° bs . This effect is more important in the Q2206A damped 
Lya system than the heating heating efficiency effect, which is why 
jstars j g higher for the "WD-high" model. Note, J^ tars for the 
"WD-high model" is smaller than for the "WD-low model" for every 
other damped Lya system in our sample. 



12 



STAR FORMATION IN DAMPED Lya SYSTEMS 



tensities derived with the C II* technique for any damped 
Lya system in our sample. This may explain why Q2206A 
is the only confirmed damped Lya system yet detected in 
emission at high rcdshift. 

4.2. J v Inferred from Direct Detection of Starlight 

M0ller et al. (2002) obtained a STIS image of the 
Q2206— 19 field with a very wide V filter in which the ef- 
fective central wavelength is 5851.5 A and FWHM=4410.3 
A. The image has an angular resolution of 0.08 arcsec. The 
field shows a diffuse galaxy with an irregular shape and a 
bright knot from which Lya emission is detected. The 
identification of the galaxy with the damped Lya system 
follows from the redshift of Lya emission, which is 250 ± 
50 km s _1 larger than the absorption-line redshift. With 
an angular separation of 9b— 0.99 arcsec, the impact pa- 
rameter of the knot with respect to the QSO sightline is 
given by 6=8.4 kpc for the WMAP (Bennett et al. 2003) 
cosmological parameters adopted here. From their STIS 
photometry of the galaxy, M0ller et al. (2002) found V = 
24.69±0.07 within a circular aperture radius of 0.3 arcsec, 
where V is the "V" magnitude corresponding to a filter 
with central wavelength A = 5851.5 A. More recently, [O 
III] 5007 emission with a rcdshift similar to that of Lya 
emission was detected from the object labeled "ps" , which 
is located within « 0.8 arcsec of the bright knot (Warren 
2004). 

To test the predictions of the C II* model we determine 
J v at the the location of the absorbing gas. Summing over 
pixels coinciding with the associated galaxy in the STIS 
image, we find 

1 1 

where we assume the physical distance between each pixel 
and the absorbing gas equals its projection on the sky. 
Here F* o is the flux density in the i th pixel, which is sepa- 
rated from the brightness centroid of the QSO by the angle 
9i and detected at frequency vq=v{1 + z) _1 , where the rest 
wavelength corresponding to v is A = 2000 A. In the pres- 
ence of dust, the sum in equation (6) would contain the 
extinction correction, exp[r|j — rj_], where Ty and t]_ are 

the 2000 A dust optical depths along the line of sight to 
pixel i and between pixel i and the absorbing gas in the 
plane of the galaxy. Independent estimates of Ty indicate 
that it varies with 9i. The depletion pattern implied by 
the ratio [Fe/Si] gas suggests rjj w 0.05 to 0.10 along the 

QSO sightline. To estimate rjj at the Lya -emitting bright 
knot, we refer to observations of Lyman Break Galaxies 
(LBGs). Shapley et al. (2003) show that E B -v for LBGs 
to decrease with increasing Lya equivalent width, W . 
Since W = 88 A, we infer Eb-v~ 0.06, and as a result 
rjjw 0.5 at the position of the knot. However, owing to the 
irregular light distribution, the main contribution to the 
sum comes from pixels closer to the QSO sightline than 
the knot; i.e., from locations where rjj << 0.5. Since t]_ 

<<0.5 in the limit 9 l — ► 0, the dust correction to equation 
(6) can be safely ignored. 

The expression for J v in cq. (6) was kindly evaluated 
for us by P. M0ller (2004), who summed F%JQ\ from an 



inner radius surrounding the QSO, 9 m i n , to an outer ra- 
dius, max . We found the result to be insensitive to the 
choice of 9 max for 9 max > 0.75 arcsec, since irregularities 
in the outer brightness distribution such as the bright knot 
contributes less than 10 % to the total J„. For this rea- 
son we let Omax = 1.5 arcsec. But the results are sensitive 
to the value of 9 m i n as it approaches the radius of the 
QSO PSF. We are confident the result is not dominated 
by PSF subtractions errors when mi „=O.36 arcsec. In 
this case we obtain a lower limit of JP hot =3.8x 10~ 18 ergs 
cm~ 2 s _1 Hz -1 sr _1 , where we have introduced the quan- 
tity JP hot to avoid confusion with the C II* technique. On 
the other hand when 9 m i n < 0.21 arcsec, the result is com- 
pletely dominated by PSF subtraction errors. Therefore, 
we obtain the "upper limit" JP hot =7.Qx 10~ 18 ergs cm~ 2 
s _1 Hz -1 sr _1 when 9 min = 0.21 arcsec. The latter re- 
sult is not a true upper limit because point sources hidden 
within the PSF could raise JP hot significantly. For exam- 
ple, if the bright knot were located within 0.1 arcsec of the 
QSO sightline, its contribution to JP hot would be a factor 
of 10 larger than that of the entire galaxy. We shall ignore 
this improbable configuration and assume a uniform light 
distribution at 9 < 9 m i n . From our solution to the transfer 
equation in WPG (see their eqs. (14) and (18)) we find 
comparable contributions to JP hot from sources on either 
side of the 9 = 9 min boundary. For this reason we con- 
clude that JP hot =(3.8^14.0)xl0- 18 ergs cm" 2 s" 1 Hz" 1 
sr _1 . 

4.3. Comparison Between C II* and Photometric 
Estimates of J„ 

Let us compare the photometric estimate of J v with the 
predictions of the C II* technique. From § 4.1 we have 
J™* =(1.7± 2 ;£)xl0- 18 ergs cm" 2 s" 1 Hz Jl sr" 1 for the 

CNM model, and J^ ir =(38±ig)x lO" 18 ergs cm" 2 s" 1 
Hz -1 sr _1 for the WNM model, where in the present sub- 
section we have substituted the quantity 11 for Jf, tars 
to avoid confusion with the photometric results (see Ta- 
ble 3). Since these intensities are evaluated at A=1500 A, 
and JP hot obtained from the photometric technique corre- 
sponds to A=2000 A, the spectral shape of the FUV con 
tinuum could be a source of additional errors. But such 
errors are negligible since the continuum is observed to be 
flat (M0ller et al. 2002), and theoretical models of star 
forming regions predict the continuum to be flat (Bruzual 
& Chariot 2003). These estimates of j£ u are compared 
with the estimates of JP hot in Table 3. Because the lat- 
ter estimates were made by equating the projected and 
physical separations between the galaxy and the absorb- 
ing gas, the results for JP hot are upper limits. As a result, 
the limits on JP hot are consistent with the values of J^ u 
predicted for the CNM, but not the WNM model. The 
discrepancies with the WNM model are larger than may 
be apparent because the error bars on J^ u correspond to 
2-a upper and lower limits set by estimates of the C II* 
column density rather than 1-er random errors obtained 
from the uncertainties in iV(H I) and to a lesser extent, 
N(C II*). 

The results in Table 3 are also illustrated in Figure 
7, which compares the photometric upper limit on JP hot 
(horizontal line) corresponding to 9 min — 0.21 arcsec with 
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Table 3 

Mean Intensities from Photometric and C II* 
Techniques 

stars d 



Technique Approximation J v " 



Notation 



Photometry 

Photometry 

C II* 
C II* 



"mill 

=0.36 arcsec 



=0.21 arcsec 
CNM 
WNM 



<3.8° 
<14.0 b 

1 7+2.7 c 

38 



j-phot 
j-phot 



j; 



+ 52 



or 

V 

tCII* 



"Inner radius used to obtain J v from eq. (6). Outer radius, ma x 
=1.5 arcsec. 

95 % confidence upper limit obtained by equating physical and 
projected separations between galaxy pixels and the QSO sightline 



c 95 % confidence errors 
J 10 



-18 



ergs cm 



Hz 



-1 



-n 



Jy predicted by the CNM model (empty triangle with 
vertical error bars). The true value of JP hot , {J^ hot )true, 
depends on the inclination angle, i, of the galaxy to the 
line of sight such that (JP hot ) true =cos 2 (i)JP hot . From Fig- 
ure 7 we see that compatibility between (JP hot )true and 
J^ n is obtained for 0° < i < 77°; i.e., for a wide range of 
plausible inclination angles. In the case of the 'BT' model 
and J bkd > 1.7(J^ kd )HM, an inclination angle, exceeding 
85° is required, which has a probability less than 0.7 %. 
Nevertheless, the uncertainties are sufficiently small for us 
to draw an important conclusion: we find a self-consistent 
model in which [C II] 158 /jm cooling is balanced by the 
heating rate supplied by the FUV radiation field of the 
galaxy associated with this damped Lya system, provided 
C II* absorption arises in CNM gas. 

5. ARE ALL DAMPED Lya SYSTEMS INTERNALLY 
HEATED BY RADIATION ? 

If the Q2206A damped Lya system is heated by local 
starlight, it is plausible to ask whether all damped Lya 
systems are heated in the same manner. This is a natu- 
ral question to pose for the damped Lya systems with 
positive detections of C II* absorption since internal heat 
sources are required in most cases (see § 3.2). Indeed WPG 
assumed the presence of internal radiation fields to deduce 
jstars £ Qr their sample of 33 damped Lya systems. While 
they ignored the background contribution to J Vl the l c (n) 
curves in Figure 3 show that background heating can be 
appreciable, especially for the "WD-high" and "BT" mod- 
els. Therefore in § 5.1 we recompute J^ tars in the presence 
of background radiation for damped Lya systems with 
positive detections. We shall also examine the role of local 
radiation for damped Lya systems with upper limits on C 
II* absorption. In § 3.2 we showed that background heat- 
ing alone could account for [C II] 158 /xm emission from 
these objects if the true value of l° bs were sufficiently low. 
In § 5.2 we shall determine whether the upper limits are 
also consistent with the presence of local radiation fields. 

5.1. Systems with Positive Detections 
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1 2 3 

J bkd /(J bkd ) H M 

solutions for Q2206A damped Lya system ver- 
sus Jl kd /(Jt kd ) H M where (Jt kd ) H M is Haardt-Madau (2003) back- 
ground radiation field used throughout this paper and J^ tars is in 
units of 10 — 18 ergs cm -2 s — 1 Hz -1 sr — 1 . Curves are dot-dashed 
for the "BT" model, dashed for the "WD-high" model, and solid 
for the "WD-low" model. Unfilled triangle labeled "J CI1 " is our 
estimate of J^ tars for "WD-low" model, where the data point is an 
average of the maximal and minimal depletion models. Error bars 
include systematic errors due to uncertainties in k, P eq , and J^ kd 
(see text). The "WD-low" solid curve was computed for minimal 
depletion. Horizontal line denoted by J phot (8 m i„=0.2l arcsec) is 
upper limit on J phot inferred from photometry of associated galaxy 
when m i„=O.21 arcsec (see text). 



We obtained J^ ars for 23 damped Lya systems with 
positive detections of C II* absorption (see Table 2). In 
each case we followed the same procedures used in § 4.1, in 
particular we added J bkd &nd J^ aTS to compute the total 
mean intensity, J u , illuminating the gas, and then con- 
sidered only solutions in which C II* absorption arises in 
CNM gas at the thermally stable density, /jcnm ■ Currently 
detected C II* absorption is unlikely to arise in WNM gas 
for most damped Lya systems with positive detections 
(see § 4). In the case of the "WD-low" CNM model the 
jstars d er i vec [ w ith the addition of background radiation 
was typically 0.1 dex less than the Jf, tars derived without 
background radiation. On the other hand, our calcula- 
tions admitted no CNM solutions for Jf / tars in the damped 
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Lya systems toward Q2241+13 and Q1036— 22: this is 
not surprising in view of the significant contribution of the 
high redshift CMB radiation to the l c {n) solutions shown 
in Figure 3b. The largest discrepancies for damped Lya 
systems that admitted CNM solutions were Q1346— 03 and 
Q0255+00 for which J^, tars with background radiation was 
0.55 dex and 0.15 dex less than without background radi- 
ation. However, these four damped Lya systems more 
closely resemble most systems with upper limits because 
l° bs is less than or equal to the peak of the predicted l c {n) 
curve. In § 5.2 we show that in systems like this, the line of 
sight likely passes through WNM gas subjected to a range 
of J^ ars including those inferred for the CNM solutions. 

Because of the low fraction of discrepant estimates of 
jstars an( j smce ^ ne discrepancies are confined to objects 
with the lowest values determined for J^ tars , we conclude 
that background radiation has little effect on the SFRs 
derived by WPG for the "WD-low" model. Specifically, 
the average Jf, tars decreased by only 0.04 dex when back- 
ground radiation was included in the models. This quan- 
tity has cosmological significance since the average J^ tars 
can be used to deduce the average SFR per unit area from 
which one can deduce the global SFR per unit comov- 
ing volume (see WGP). The same qualitative arguments 
hold for the "WD-high" and "BT" models. The average 
jstars decreased by 0.15 dex when background radiation 
was included in the "BT" model and by 0.3 dex for the 
"WD-high" model. 

5.2. Systems With Upper Limits 

The nature of damped Lya systems with upper limits 
on C II* absorption is determined by the ratio of the true 
value of l° bs to the peak value of l c {n) predicted for back- 
ground heating alone, l v c eak . This point is illustrated in 
Figure 8, which shows l c {n) curves predicted by the "WD- 
low" model for the Q2348-14 and Q2237-06 damped Lya 
systems. These objects were selected to represent the high- 
est and lowest values of the ratio l° bs /lP eak in Figure 4. We 
computed l c (n) curves for the case J u =J bkd . This results 
in the lowest l c {n) curves in Figure 8 (shown as red dotted 
curves), which correspond to background heating alone. 
We then let J^,—J bkd +J^ tars to account for the presence 
of local radiation and computed l c (n) for the same grid of 
SFRs per unit area, -4.0 < logioV'* < 0.0 M yr" 1 kpc~ 2 , 
used to find J° tars for the Q2206A damped Lya system. 
This corresponds to mean intensities in the range 5x 10~ 21 
< J? ars < 5xl0~ 17 ergs cnT 2 s" 1 Hz" 1 sr" 1 at A = 1500 
A (see WPG and § 5.1). By comparison, J bkd = 2.7xl0- 20 
and 2.0xl0~ 20 ergs cm~ 2 s _1 Hz -1 sr _1 at the redshifts of 
the Q2348-14 and Q2237-06 damped Lya systems. The 
l c {n) solutions including J^, tars are shown as solid curves 
in the Figure 8. 

First, we consider the results shown in Figure 8a for 
Q2348-14. In this case the upper limit on l° bs (10^ 26 9 
ergs s _1 H _1 ) is significantly higher than iv eak (10~ 27,4 
ergs s _1 H _1 ). If future spectroscopy shows the true value 
of l° bs to remain above l p c eak , then l° bs could intersect 
WNM or CNM thermally stable solutions corresponding 
to a wide range of J** ars : the (n,l c ) coordinates of these 
stable solutions are respectively shown as the steeply ris- 
ing dot-dashed and solid curves in Figure 8a. The crosses 
in Figure 8a depict upper limits on n set by the 



C II/C I constraints for each l c (n) curve. Because these 
density limits are less than ticnm, they appear to rule 
out the CNM solutions. However, WGP showed the lim- 
its to be sensitive to local input parameters such as X- 
ray and cosmic-ray ionization rates. For these reasons we 
make the conservative estimate that the CNM models are 
not excluded for upper limits on n exceeding sa 0.2ncNM 
Therefore, if l obs > P c eak , the CNM and WNM models are 
in principle equally plausible for this damped Lya sys- 
tem. But, because we have found direct evidence against 
the WNM hypothesis for other damped Lya systems with 
similar ratios of l obs to V c eak (see Howk et al. 2004 a,b) 
we believe the CNM hypothesis to be more likely in these 
cases. 

On the other hand, the CNM solutions for Q2348-14 
are implausible if l obs < l peak . The presence of background 
radiation causes a significant increase in l c (n) for low val- 
ues of Jf, tars . The result shown in Figure 8a is that no 
thermally stable CNM solutions are possible for l° bs less 
than the intersection of the background heated l c {n) and 
the locus of thermally stable CNM solutions; i.e., when 
l obs < 10~ 27 ' 5 ergs s" 1 H" 1 . By contrast, the WNM 
solutions, thermally stable or otherwise, are consistent 
with the condition l° bs < V^ ak for the entire grid of J* tars . 
The low densities, n < 1 cm~ 3 , predicted for models with 
large values of J^ tars argue in favor of moderate values 
of Jl tars (< 10~ 18 ergs cm~ 2 s _1 Hz -1 sr _1 ). Otherwise 
the linear dimensions of the absorbing gas clouds would 
become implausibly large. It is important to note that 
X-ray heating dominates grain photoelectric heating in all 
the WNM solutions, and that the Xray background domi- 
nates local X-ray sources for Jf, tars < 5xl0~ 19 ergs cm~ 2 
s- 1 Hz" 1 sr- 1 . ' 

To investigate the effect of background radiation on these 
conclusions, we recomputed the l c (n) curves with the same 
procedures, but excluded all radiation backgrounds. The 
results are shown in Figure 8b. Without background heat- 
ing no lower bound on thermally stable CNM solutions 
exists; i.e., the thermally stable CNM locus can extend 
to arbitrarily low values of l c {n). Therefore, when the ef- 
fects of background radiation are omitted, CNM solutions 
with l obs < F c eak can exist in principle. But when back- 
ground radiation is present and l obs < V c eak , the C II/C I 
constraint results in upper limits on n < 0.03ncNM (Fig- 
ure 8a). Because this discrepancy is larger than allowed 
by model uncertainties, we conclude that thermally stable 
CNM solutions are ruled out when l obs < V c eak . 

Next we consider the results for Q2237— 06 in Figure 8c. 
In this case the existing upper limit on l° bs is less than l"P eak 
and as a result all of the CNM solutions are ruled out. Be- 
cause of the large redshift, z=4.080, the CMB dominates 
the behavior of l c {n) for low values of Jf, tars . This is illus- 
trated in Figure 8d, which plots the l c (n) curves that result 
in the absence of background radiation. Comparison with 
Figure 8c reveals two features caused by background radi- 
ation: a low-density asymptote at l c (n) — 10^ 27 - 6 ergs s -1 
H _1 due to CMB radiation, and an increased peak value, 
lpeak = ^Q-27.3 er g S g -i jj-i^ j ue ^ Q ^ e combination of 

CMB and X-ray background radiation. Again, Figure 8d 
shows that without background radiation, thermally sta- 
ble CNM solutions would be possible. Note, the absence 
of density limits is due to the lack of a C II/C I constraint 
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Fig. 8. — Solid curves are equilibrium £ c (n) solutions for gas subjected to internal heating rates generated by a grid of SFRs per unit area, 
logio^*=-4.0,-3.5,....,O.OM Q yr" 1 kpc~ 2 . Same form of FUV, X-ray, and cosmic-ray inputs as in Figure 6. (a) and (b) show results for the 
Q2348— 14 damped Lya system in which background radiation is included in (a) and excluded in (b). Lowest dotted curve in (a) is the result 
for background radiation alone. The horizontal line with arrow corresponds to 2— a upper limit on l° bs . The steeply rising dot-dashed and 
solid curves are thermally stable WNM and CNM solutions at n=»wNM an< i n CNM- Crosses are upper limits on n set by C II/C I constraint 
for each (c) and (d) are analogous results for Q2237— 06 damped Lyc* system, (see § 4.2 for discussion) 



for this damped Lya system. 

Therefore, since l obs jl v ^ ak < 10 02 for most cases in Fig- 
ures 4a and 4b, we predict the sightlines through most 
damped Lya systems with upper limits on C II* absorp- 
tion to pass through WNM gas alone. The same con- 
clusion holds for the damped Lya systems with positive 
detections for which l° bs jl v c eak are also less than 10 02 (i.e., 
the damped Lya systems toward Q0127-00, Q0255+00, 
Q1036-22, Q1346-03, and Q2241+13). In the few cases 
of upper limits with higher ratios of l obs to lP eak , the sight- 
lines could encounter either WNM or CNM gas. But the 
similar l° bs jl v c eak ratios for damped Lya systems with pos- 
itive detections, objects for which there is independent ev- 
idence favoring the CNM hypothesis (cf. WGP and Howk 
et al. 2004), suggests that in these cases the gas is CNM. 

6. THE IMPACT OF IONIZED GAS ON THE CII* ANALYSIS 

The thermal equilibria presented above were computed 
by assuming C II* absorption arises in neutral gas. But, 
C II* absorption could also arise in ionized gas. Since 
a significant fraction of the C II* absorption detected in 
the halo and thick gaseous disk of the Galaxy arises in a 
warm ionized medium (WIM; see Lehner, Wakker, & Sav- 
age 2004), the same may be true in damped Lya systems. 



Indeed, the presence of Al III, C IV, and Si IV absorption 
in most damped Lya systems indicates that the QSO 
sightlines typically pass through ionized as well as neutral 
gas (Wolfe & Prochaska 2000a,b). 

However, there are several reasons why C II* absorp- 
tion in damped Lya systems is unlikely to arise in ionized 
gas. First, from the constraint on the optical depth ratio, 
r(Si II* 1264.7) /t(C II* 1335.7) <1(T 2 - 5 , Howk, Wolfe, & 
Prochaska (2004a) deduced an upper limit of T < 800 K 
for the z=4.22 damped Lya system toward PSS 1443+27. 
C II* absorption cannot arise in ionized gas in this damped 
Lya system because such gas typically has temperatures 
exceeding 10 4 K. Because there are indications of similar 
limits in at least two other damped Lya systems, this 
temperature limit may be a generic property, which would 
rule out ionized gas as the site of C II* absorption in most 
damped Lya systems. Second, consider the highly ion- 
ized gas, which gives rise to C IV and Si IV absorption. 
In a systematic study of 32 damped Lya systems, Wolfe 
& Prochaska (2000a, b) found statistically significant dif- 
ferences between the velocity structure of the low ions and 
the high ions. Because the C II* velocity profiles are, with 
rare exceptions, indistinguishable from the low-ion reso- 
nance line profiles (WPG), C II* absorption is unlikely to 
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originate in highly ionized gas. This has interesting impli- 
cations. The presence of narrow C IV absorption com 
ponents with FWHM < 30 km s" 1 implies T < 2xl0 5 
K, which indicates CIV is photoionized in damped Lya 
systems. To summarize, although hot gas may be present, 
comparison between the velocity structure of the C IV and 
C II* absorption profiles suggests the two ions do not co- 
exist in the same gas, and the narrow C IV line widths 
indicate C IV is not directly tracing hot gas nor its inter- 
action with cool gas. 

On the other hand, Lehner, Wakker, & Savage (2004) 
present convincing evidence that a significant fraction of 
C II* absorption lines in the ISM arises in the WIM. Since 
the WIM of the ISM of the Galaxy is also the site of Al III 
absorption (Savage & Edgar 1990), the implication is that 
the gas giving rise to Al III absorption is also the site of C 
II* absorption. Because Al III absorption is a generic fea- 
ture in the spectra of damped Lya systems, and because 
of the similarity between the centroids of the multiple ve- 
locity components in the Al III and low-ion profiles (Wolfe 
& Prochaska 2000a), the close resemblance between the C 
II* and low-ion profiles suggests that C II* absorption in 
damped Lya systems also arises in a WIM. However, this 
is unlikely. First, in some instances the ratio of Al III 
to low-ion column densities varies strongly with velocity 
(Wolfe & Prochaska 2000), indicating that Al III may not 
be an accurate tracer of C II* absorption in damped Lya 
systems. Second, even if Al III does trace the low ions, 
Al III need not trace WIM gas in damped Lya systems, 
since there are other means of producing Al III than stel- 
lar photons. At high rcdshifts a primarily neutral gas can 
contain significantly higher column densities of Al III than 
at z—Q; i.e., at high z the presence of Al III need not be a 
signature of a WIM. This is due to the predicted increase 
with z of the soft X-ray background intensity (Haardt & 
Madau 2003): at z > 2 the background intensity at 0.5 keV 
(shown in Figure 1) should be more than 30 times higher 
than at z = 0. As a result, more radiation is available 
for photoionizing Al II (IP=18.8 cV) to Al III (IP=28.4 
eV) than at z=0. Prochaska et al. (2002) find empirical 
evidence for this in the z=2.62 damped Lya system to- 
ward Q1759+75, where comparison between intermediate 
and low-ion velocity profiles shows significant Al III absorp 
tion at the position of the velocity component associated 
with neutral gas. 

These points are illustrated in Figure 9, which shows 
the results of a CLOUDY (Ferland 2001) photoionization 
calculation in which we simulated the damped Lya sys- 
tem toward Q2206A with a slab with AT(H I) = 10 20 - 5 
cm~ 2 , n=0.1 cm~ 3 , and [Fe/H]=— 1. The slab is sub- 
jected to the same Haardt-Madau backgrounds (2003) as 
shown in Figure 1. While local inputs were neglected, 
these will have negligible effect on the Al II/ Al III ratio 
or the (CI+CII)/C ratio. The x axis is the ionization pa- 
rameter, the blue curve shows the [Al III/A1 II] ratio, the 
red curve shows the linear neutral gas fraction, 1—x, and 
the orange curve shows logi (CII/C). Since we have not 
measured Al II 1670 (which is undoubtedly saturated) we 
assume [Al II/Fe II] =0 which is found for the few damped 
Lya systems in which Al II 1670 is unsaturated. In that 
case the observed ratio, [Al II/ Al III] = 0.58 dex. The 
dotted vertical lines show that value ± 0.1 dex. The cor- 



responding neutral gas fraction, 1—x, runs from 0.86 to 
0.95. Therefore, one can accommodate a fairly large Al 
III column density in a mainly neutral gas, which pro- 
duces significant amounts of C II. Furthermore, Figure 9 
reveals a subtle but important point on the multi-phase 
nature of the damped Lya system. If one assumes a single 
phase, then the gas is w 10% ionized and cannot be con- 
sidered a cold neutral medium. This result contradicts our 
assessment of the C II* profile and we are driven to a two- 
phase solution (WNM+CNM) in order to accommodate 
the A1III and CII* measurements. Assuming compara- 
ble HI column densities for the WNM and CNM phases, 
we find a solution with logt/ w —4.5 for the CNM and 
log U w —3 for the WNM, which gives an average Al II/ Al 
III ratio consistent with our observations. Adopting J VH 
w 10~ 21,5 for the ionizing background radiation, we find 
the implied density values are consistent with the model 
presented in Figure 6, e.g., logio^CNM ~ 0.2cm -3 . Note 
that this relatively low density for the CNM is due to the 
high metallicity of the system. In summation, the pres- 
ence of Al III does not require a WIM and in fact the 
observed Al II/ Al III ratio is consistent with one's expec- 
tation based on the two-phase models derived from the 
CII* observations. 

Further evidence against the ionized gas hypothesis for 
C II* absorption can be found in Figure 10, which plots 
l° bs versus N(H I) for the sample of damped Lya systems 
in Table 2. The data show no evidence for correlations 
between l° bs and N(R I): the Kendall tau test reveals the 
probability for the null hypothesis of no correlations to be 
65 %. By contrast, in a recent study of C II* absorption in 
the ISM, Lehner, Wakker, & Savage (2004) find a system- 
atic increase in l° bs with decreasing N(H I), which they 
reasonably attribute to an increase in ionization level in 
systems with low iV(H I). To understand this we note that 
in the absence of radiative excitation, l° bs =nAcn(l— a;) -1 , 
where nAcn is the 158 /Kin cooling rate per particle . In 
the low-density, high-ionization limit applicable here, l° bs 
oc J VH , where J VH is the radiation mean intensity at the 
Lyman limit. As a result, l° bs will increase with increasing 
J VH , which causes N(H I) to decrease for a 
fixed total column density. By contrast, we find no con- 
vincing evidence for such an anti-correlation at JV(H I) > 
2xl0 20 cm -2 , which is additional evidence supporting the 
hypothesis that C II* absorption in damped Lya systems 
arises in mostly neutral gas. 

Finally, our arguments against background radiation as 
the heat source source for damped Lya systems with pos- 
itive detections of l° bs depended on the upper limits on 
the total density, n, inferred from photoionization equilib- 
rium between C II and C I. The C II/C I ratios and the 
resultant density limits were computed for gas that was 
mainly neutral, and in the WNM phase for most cases. 
While ionized gas may not be the site of C II* absorption, 
it may still be present and could contribute to the total 
N(C II). Taking the results discussed above for Q2206A, 
in which x < 0.14, and assuming N(C V)—0 and N(C 
II)=(C/H)iV(H)v7iM for the WIM (where N(R) is the to- 
tal column density of H), then the total C II column den- 
sity, N(C II)=(C/R)(N{H) wim +N(R) WN mM1/0.86) 
x(C/H)AT(H I) = 1.16(C/H)V(H I). Since the results are 
the same for N(C I), the change in the C II/C I ratio will 



Wolfe, Prochaska, & Gawiser 



17 



> 



1.4 

1.2 

1.0 

0.8 

0.6 

0.4 

0.2 
0.0 



i — i — i — i — i — i — r — i — i — I — i — i — i — i — i — i — i — v. — i — I — i — i — i — i — r — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — i — i— 




log[N(AlII)/N(AHII)] 
1-x 

(CII+CI)/C 




-6 -5 -4 -3 -2 -1 

logU 

Fig. 9. — Various ionization ratios versus ionization parameter resulting from simulation of photoionization equilibrium for the Q2206A 
damped Lya system. We let ionizing background radiation be incident on a plane parallel slab with N(H I)= 10 20 ' 5 cm -2 , n=0.1 cm -3 , 
and [Fe/H] = —1. Blue curve is [A1II/A1III], red curve is linear neutral gas fraction, 1— x, and orange curve is (CII+CI)/C. 



be at most 16%. This is insufficient to alter the conclusions 
given in § 3.1. 

To summarize, by contrast with the ISM, the evidence 
accumulated so far suggests that the WIM does not give 
rise to significant C II* absorption in damped Lya sys- 
tems. While an investigation of the difference in physical 
conditions is beyond the scope of this paper, we suggest 
two possibilities: (1) an extensive WIM may not be present 
in damped Lya systems owing to a low escape fraction 
of ionizing radiation from H II regions embedded in the 
neutral gas. This is consistent with the presence of Al III, 
the main indicator of the WIM in the ISM, since Al III 
can arise in neutral gas in damped Lya systems due to 
the increased strength of the X-ray background at large 
redshifts. (2) If an extensive WIM does exist in damped 
Lya systems, C II* absorption is suppressed because of 
the low value of (C/H) and possibly because the mean elec- 
tron density is lower than in the neutral ISM (see eq. 4). 
Lower electron densities could result from gas pressures, 
which are lower than in the CNM/WNM interface. 

7. CONCLUDING REMARKS 
Our main conclusions are: 

(1) Most damped Lya systems with positive detections 
of C II* A 1335.7 absorption are not heated by external 
background radiation, but rather require an internal en- 
ergy source. 

(2) The internal energy source for one such object, the 
Q2206A damped Lya system, appears to be FUV radia- 
tion emitted by stars in the associated galaxy identified by 
M0ller et al. (2002). In this case, C II* absorption cannot 



arise in WNM gas, but rather arises in CNM gas. 

(3) Independent arguments, in particular the strong cor- 
relation between l° ba and dust-to-gas ratio, suggest C II* 
absorption currently detected in most damped Lya sys- 
tems also arises in CNM gas that is similarly heated. 

(4) Damped Lya systems with upper limits on C II* 
absorption also require internal energy sources if the true 
values of l° bs exceed the peak values of l c (n) predicted for 
heating by background radiation alone, l^ eak . While the 
absorbing gas could in principle be either CNM or WNM, 
the resemblance of the ratio l° bs /lP eak to most systems 
with positive detections suggests the gas is CNM. 

(5) For most damped Lya systems with upper limits 
the true values of l° bs are likely to be less than l v c eak . In 
this case the gas could be heated by background radiation 
alone. But the data are also consistent with the added 
presence of internal energy sources. In either case, ab- 
sorption is likely to arise in WNM gas. The same conclu- 
sion holds for the few damped Lya systems with positive 
detections and in which l° bs is also less than l v c eak . 

(6) In the case of radiative heating, the radiation inten- 
sities deduced for the positive detections and allowed by 
the upper limits are about the same; i.e. 10~ 19 < Jf ars 
< 10~ 18 ergs cnr 2 s _1 Hz -1 sr _1 . If the radiation is 
emitted from a uniform distribution of stars in a disk, the 
range of SFRs per unit area is given by 10 -3 < -0* < 10 -2 
Mq yr _1 kpc -2 (see WPG). Therefore, all damped Lya 
systems may be drawn from a single parent population of 
two-phase gas layers for which the strength of C II* ab- 
sorption depends only on the phase encountered by the 
line of sight through the damped Lya system. In the case 
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Fig. 10. — l° ba versus JV(H I). Solid red, green, and blue data 
points are positive detections, 95 % confidence upper limits, and 95 
% confidence lower limits for damped Lya systems in Table 2. 



of strong C II* absorption, the gas includes a CNM with 
T ~ 100 K and n ~ 10 cm -3 , and for weak absorption the 
gas is WNM alone with T~ 8000 K and n ~ 10" 1 cm" 3 . 
Therefore, detectability of C II* absorption appears to de- 
pend only on the density of the gas encountered by the 
line of sight. 

(7) By contrast with the ISM, the evidence accumulated 
so far suggests that C II* absorption in damped Lya sys- 
tems does not arise in a WIM. 

The arguments against background radiation as the dom- 
inant source of heat input for damped Lya systems with 
positive detections of C II* absorption are robust as they 
are insensitive to the properties of the dust responsible for 
grain photoelectric heating. Because of the low values of 
J^, kd , the large lower limits on C II/C I imply small upper 
limits on n, typically n < 0.03 cm~ 3 (sec equation 5). In 
that case the gas is restricted to the WNM phase where 
the predicted 158 /mi emission rate, l c (n), is not tied to the 
grain photoelectric heating rate. Instead, l c (n) depends on 
n e , which is well determined for the input quantities n and 
J^, kd , and the fine-structure collision strength, which is ac- 
curately determined. We also find that l c (ri) in the WNM 
is insensitive to variations in background intensity: when 
we increased J^, kd by a factor of 10, l c (n) in the WNM in- 



creased by only 0.2 dex. However, to relate temperature to 
density we assumed thermal equilibrium, which may not 
be applicable owing to the relatively long cooling times of 
the WNM (cf. Wolfire et al. 2003). But it is difficult to 
envisage any scenario in which 158 /jm emission from gas 
at such low densities depends on grain properties. 

By contrast, the conclusions concerning the heat source 
for the Q2206A damped Lya system are sensitive to the 
nature of the dust responsible for grain photoelectric heat- 
ing, because it is so metal rich. Specifically, J^ tars for the 
"BT" model is between factors of 3 to 5 lower than for the 
other models. The differences are related to differences in 
heating efficiencies, and differences in the predicted values 
of ncNM- Despite this, the agreement between our esti- 
mates of JP hot and Jf ars (illustrated in Figure 7) indi- 
cates that FUV starlight emitted by the galaxy associated 
with this DLA is incident on the grains associated with 
the gas and heats the gas by the grain photoelectric mech- 
anism. This further argues against mechanical heating 
mechanisms such as turbulent dissipation (see discussion 
by Wolfire et al. 2003). We emphasize that the agreement 
holds only for the CNM model, since the WNM predictions 
for all three grain models indicates J^* ars exceeds the 95 
% confidence upper limit on JP hot . This result is in agree- 
ment with the conclusions reached by WPG and WGP 
who found that the bolometric background intensity pre- 
dicted at z = if the WNM model applies to all damped 
Lya systems violated observational constraints. We also 
emphasize that the J^ tars predicted for this damped Lya 
system are the highest in our sample, which may help to 
explain why Q2206A is the only confirmed damped Lya 
system detected in emission at high redshift. 

Finally, we consider the SFR of the Q2206A damped 
Lya system. Summing over the same pixels used to de- 
termine Jf, tars in equation (6), we can estimate the total 
FUV flux density of the associated galaxy. Assuming again 
that Bmax = 1.5 arcsec, we find that U=23.0 if f? min =0.36 
arsec and F=22.9 if 9 min = 0.21 arcsec; i.e., the results 
are insensitive to the value of m i n . We also find they 
are insensitive to the value of 9 ma x for 6 max > 1.0 arc- 
sec. Assuming U=23 we find that the luminosity per unit 
bandwidth, i^=(2.0±0.50) x 10 29 ergs s _1 Hz -1 where we 
have adopted conservative 25 % errors for the photome- 
try. As a result SFR = (26±6.5)M© yr _1 . Because we 
have not corrected for the effects of extinction, this is a 
lower limit. Following the discussion in § 4.2 in which 
we used the equivalent width of Lya emission to esti- 
mate Eb-v we find that SFR < 50 M Q yr^ 1 . Comparison 
with LBGs, which are drawn from a Schechter luminosity 
function characterized by (SFR)*=60 M Q yr _1 and a = 
1.6 (Steidel et al. 1999; Shapley et al. 2003) shows the 
Q2206A galaxy to be an LBG with a typical SFR. On 
the other hand its irregular and extended morphology is 
atypical for LBGs and may be atypical for galaxies associ- 
ated with damped Lya systems. For this reason it is not 
obvious that these results can be extrapolated to other 
damped Lya systems with detected C II* absorption by 
simply scaling SFR with J^ tars . Indeed, applying the sim- 
ple uniform disk model of WPG and WGP to Q2206A, we 
find tp* — 0.03 M© yr _1 kpc~ 2 rather than the observed 
value of > 1 M© yr _1 kpc -2 . The latter value of ip* 
is consistent with the derived value of JE hot because of 
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the small solid angle subtended by the FUV continuum 
sources from the location of the absorbing gas. Of course 
the SFR per unit H I area would be considerably lower 
if the star forming region were embedded in a larger H 
I envelope. Interestingly, at least 50 % of the remaining 
damped Lya systems with detected C II* absorption are 
predicted to have J° tars within a factor of 3 of the Jf ars 
inferred for the Q2206A damped Lya system. It is not 
clear at this time whether or not their SFRs are within a 
factor of 3 of the SFR derived for the Q2206A galaxy. We 
are pursuing a program with the ACS aboard the HST to 
investigate this question. 
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